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The color neutral two-flavor superconducting (2SC) phase of cold and dense quark matter is 
studied in the presence of constant magnetic fields and at moderate baryon densities. In the 
first part of the paper, a two-flavor effective Nambu-Jona-Lasinio (NJL) model consisting 
of a chiral symmetry breaking (xSB) mass gap as, a color superconducting (CSC) mass 
gap As and a color chemical potential /ig is introduced in the presence of a rotated C/(l) 
magnetic field B. To study the phenomenon of magnetic catalysis in the presence of strong 
. magnetic fields, the gap equations corresponding to gb and As, as well as /zg are solved in 

fvj \ the lowest Landau level (LLL) approximation. In the second part of the paper, a detailed 

numerical analysis is performed to explore the effect of any arbitrary magnetic field on the 
above mass gaps and the color chemical potential. The structure of the xSB and CSC phases 
is also presented in the ^^ — eB plane, and the effect of /xg on the phase structure of the 
model is explored. As it turns out, whereas the transition from the xSB to CSC phase is of 
first order, nonvanishing /ig affects essentially the second order phase transition from CSC 

i^ ' to the normal phase. 

JiLl' PACS numbers: ll.30.Qc, 12.39.-x, 12.38.-t, 12.38.Aw 

(N 

> ■ I. INTRODUCTION 

(N 

<^ • Dense baryonic matter at low temperature and asymptotically large chemical potential is known to 

be a color superconductor [1]. This can be shown in the framework of perturbative Quantum Chro- 

(^ • modynamics (pQCD). To explore the color superconducting phase at moderate chemical potential, 

^—^ ■ however, it is necessary to use effective models, such as the well-known NJL model with four-fermion 

interaction [2\. Using an appropriate NJL type model, one can show that at baryon densities /Xc — 350 

vj MeV, i.e. only several times larger than the density of nuclear matter, the two-flavor color super- 

s' ■ conducting (2SC) phase might be present [3|, |j] (see [5(] for recent reviews on color superconductivity 

in dense quark matter). Different astrophysical processes might therefore be influenced by the color 
superconductivity that is supposed to exist inside the compact stars. In [(|, the competition between 
the chiral symmetry breaking and the color symmetry breaking condensates is investigated in the 
framework of a two-flavor color neutral NJL type model, including meson and diquark condensates, 
(To and Ag. Imposing the color neutrality condition, it is found that in the 2SC phase at /x > /x^ = 342 
MeV, the color chemical potential /xg acquires rather small values of about 10 MeV.^ Here, /ic is the 
critical chemical potential. The diquark mass gap is numerically computed to be Aq — 100 MeV. It is 
also shown that the appearance of a coexistence regime (mixed phase) depends directly on the relative 
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strength of the meson and diquark couphng constants Gs and Gd- This is also indicated in 
where it is stated that neglecting the quark masses and choosing Gd < Gs, no mixed phase appears 
at /i > fie- The xSB and CSC phases can therefore be studied separately under these conditions. 

In the present paper, we study the mesons and diquarks in the color neutral 2SC phase of cold 
and dense quark matter in the presence of constant magnetic fields. The aim is to study the effect of 
the magnetic field on the formation of chiral as well as diquark condensates, the dependence of mass 
gaps on the chemical potential // and the external magnetic field, the phase diagram // vs. B, and the 
effect of nonvanishing color chemical potential on the type of phase transitions for different fx and B 
at zero temperature T? 

The study of quark matter in the presence of constant magnetic field is relevant for the astrophysics 
of compact stars: Strong magnetic fields exist on the surface of compact stars. For neutron stars the 
magnetic fields B < 10^^ GauB, whereas for magnetars, they can be as large as i? ~ 10^^ Gaufi [12l |. 
n the interiors of compact stars, the magnetic field can be even several orders of magnitude larger 
131]. On the other hand, it is believed that the superdense interior of compact stars may be composed 
of electric and color neutral quark matter in the color superconducting phase. To test the predictions 
of astrophysical signatures of color superconductivity, a better understanding of the role of magnetic 
fields on the CSC phase is important. The study of superconducting phase in the presence of external 
magnetic fields is also relevant for the physics of heavy ion collisions: According to [ij, ll5|, in off- 
central collisions, heavy ions possess a very large angular momentum and very strong magnetic fields 
can be created. In [iq], it is shown that the magnetic field presently created at RHIC is at most 
eB ~ 1.3 m^ ~ 4.3 x 10^^ GauB, and the estimated value of the magnetic field strength for the LHC 
energy amounts to 15 m^ ~ 5 x 10^^ Gauf^.'^ Recently, the question of accessibility of the 2SC phase 



in the future heavy ion collision experiments is investigated in IT]]. Here, the authors do not consider 
the effect of the before mentioned magnetic fields. It would be therefore important to study the effect 
of external magnetic fields on the formation of 2SC diquark condensates, as well as the corresponding 
phase structure in the presence of external magnetic fields. As for the results presented in this paper, 
they may be relevant only for the physics of the heavy ion collisions, because in contrary to the electric 
and color neutrality requirement of the superdense core of the compact stars, only the color neutrality 
condition is considered in this paper. 

The effect of constant magnetic field on the formation of diquark condensates has been investigated 
by several authors. In [ig . Il9t |. it is shown that there is a linear combination of photon and a gluon 
that remains massless. The resulting "rotated" external magnetic field can therefore penetrate the 
color superconducting region without bein g af fected by the Meissner effect. This has consequences for 



the structure of compact star cores. In [20|, |2l[, the formation of magnetic color-flavor locked (MCFL) 



phase, as well as the transition to the paramagnetic-CFL (PCFL) phase are studied. In [22l. |23|]. it is 
shown that for small magnetic fields, the CFL mass gap as well as the corresponding magnetization 
exhibit small oscillations, the van Alfven-de Haas (vAdH) oscillations. This effect, which is well-known 
from condensed matter physics, is predicted by Landau and observed experimentally by van Alfven 



and de Haas (see 2j] for an investigation of this effect in cold dense quark matter in a homogeneous 



^ The effect of finite temperature will be presented elsewhere [11 
^ Here, the pion mass, m^ — 140 MeV. 



magnetic field). The transport properties of 2SC phase is investigated recently in 25l |. 



Recently, in [26|], the formation of chiral and diquark condensates as well as their competition in 
the 2SC phase at zero temperature and moderate densities are studied using the same NJL type model 
as in the present paper. It is shown that for vanishing magnetic field, a mixed broken phase can be 
found where both chiral and superconducting gaps are non-zero. For eB = 0.05 GeV^ (corresponding 
to i? ~ 8.5 X 10^*^ Gaufi) and moderate diquark-to-chiral coupling ratios Gd/Gs, the chiral and 
superconducting transitions become weaker. For large Gd/Gs, strong magnetic fields disrupt the 
mixed broken phase region and a first order phase transition is found between the xSB and the CSC 
phase for eB = 0.05 GeV^. In contrast to [26|], our results include a detailed analytical and numerical 
survey on the effect of external magnetic field and color chemical potential on cold and dense as well 
as color neutral quark matter in the presence of external magnetic fields. 

The organization of this paper is as follows: In Sec. II, starting from an extended Lagrangian density 
of a gauged NJL model containing two flavors, and following the method presented in jlSl . [lOj, we 
introduce the rotated magnetic field eB and determine the Lagrangian density containing the xSB and 
CSC mass gaps, as and A^, respectively. In Sec. Ill, the one-loop effective action and thermodynamic 
potential of the model are determined at zero temperature and finite quark chemical potential. In Sec. 
IV, assuming very strong magnetic fields, we solve analytically the gap equations corresponding to as 
and Ab, as well as the color chemical potential /ig in an appropriate LLL approximation. The xSB 
and the CSC phases are studied, in IV. A and IV. B, separately. This is possible because of our specific 
choice of free parameters, the quark mass ttiq and the meson and diquark couplings Gs and Gd- In 
the xSB phase, characterized by us / and A^ = ^g = 0, the magnetic field enhances the bound 



state formation. This is because of the phenomenon of magnetic catalysis [30|, |3l|] studied intensively 
in the past few years. ^ In the CSC phase, characterized by as = 0, A^ 7^ and /ig / 0, we determine 
analytically the fi and eB dependence of A^ and /ig in the regime of LLL dominance. In Sec. V, 
a numerical analysis is performed to study the eB dependence of the xSB and CSC mass gaps at 
// = 250 MeV (in the xSB regime) and // = 460 MeV (in the CSC regime). For small values of eB, we 
observe vAdH oscillations in the mass gaps as well as the corresponding magnetizations, as expected. 
These are also observed in [23, |23| for three-flavor CFL phase at // = 500 MeV. At eB ~ 0.45 - 0.5 
GeV^, the oscillations end up in a "linear regime". Comparing eventually our numerical results for 
eB > 0.45 GeV^ with the analytical results arising in Sec. IV for strong magnetic fields in the LLL 
approximation, we conclude that this approximation is only reliable in the above linear regime. The 
/U-dependence of the mass gaps and the color chemical potential is also discussed for various eB. 
Our results for vanishing eB coincide with the results in [6|. We also present the phase structure of 
xSB and CSC phases in a ^c — ei? plane. In particular, we are interested on the effect of the color 
chemical potential /ig on the phase structure of the model. As it turns out, for /ig = 0, a first order 
phase transition exists between the xSB and the CSC phase in the regime fic — 350 — 450 MeV and 
eB £ [0,0,7] GeV^, whereas the transition from the CSC to the normal phase is of second order and 
occurs at fic ~ 750 — 800 MeV. For //g 7^ 0, however, whereas the phase transition between the xSB 
and the CSC phase is still of first order, the second order phase transition between the CSC and the 



"* See [22] for the application of magnetic catalysis in cosmology, [33] for its application in condensed matter physics, 
and [34| - [36l . l4lj for its applications in particle physics. 



normal phase goes over into a first order phase transition between the CSC and the normal phase at 
fic — 755 MeV and eB ~ 0.13 GeV^. Note that the first order nature of the transition between the 



xSB and CSC phases was expected from 26|], where the type of phase transition between these two 
phases is studied for a fixed eB = 0.05 GeV^ and various Gd/Gs ratios. Our results confirm the 



findings in 26| for a wide range of eB G {0, 0.7} GeV^ and fixed value of Gd/Gs = 0.75. Section VI 



is devoted to a summary of our results and concluding remarks. 

II. TWO FLAVOR 2SC MODEL AT T = 0, AND /i, B 7^ 

Let us start with the fermionic part of the extended Lagrangian density of a gauged NJL model^ 

Cf = ^{x)[ijf'{d^ - ieQA^ - igT'^Gl) - mo + /i7V(^) 

+Gs[{i^{x)^{x)f + (V^(x)i75fV'(x))2] + GD[(iV^^(x)e/ee75V'(^))(iV^(a:)£/ec75V'^(x))]. (II.l) 

Here, if) = Gtp and ip = ipC are charge-conjugate spinors, and G = ^7^7*^ is charge-conjugation 
matrix, r = {ti,T2.,t^) are Pauli matrices. Moreover, (ej?)"'' = {ecY^^ and {£f)ij are antisymmetric 
matrices in color and flavor spaces, respectively. For a theory with two quark flavors, i,j = (1,2) = 
{u, d), and three color degrees of freedom a,b = (1, 2, 3) = (r, g, b). We assume that both quarks have 
the same (bare) mass ruu = ma = mo.^ Further, /t is defined by /i = ^ + -v/S/igAg, where n is the 
quark chemical potential and is responsible for the nonzero baryonic density of quark matter, and 
;Ug is inserted by hand to impose the color neutrality after the process of dynamical color symmetry 
breaking. Here, T^ = -^ with Ag = -7=diag(l, 1, — 2) the 8 Gell-Mann A-matrix. The scalar and 
diquark couplings are denoted by Gs and Gd, respectively. Furthermore, Q = Qj Ic with Qj = 
diag (2/3, — 1/3) is the fermionic charge matrix coupled to f/(l) gauge field A^^. The same setup 
without the coupling to A^j, and G^ is also used in J6J. To determine the effective action of the above 
model, we introduce first the bosonized Lagrangian density 

Cf = V^(a;)[i7^(a^ - ieQA^ - igT^G^^) + fi^^Mx) - ^(x)(m + ij^T-n)ij{x) 

1 —r 1 — 0-2 + vf^ A^A*^ 
--A*'iiij ix)efeh5i^{x)) - -A\i^P{x)efelro^'^{x)) - ^^- ^^— , (n.2) 

with m = mo + o", that includes the auxiliary mesonic fields 

a = -2Gs(?V'), 7f = -2G5(M^T^), (11.3) 

and diquarks 

A3 = -2GD{iii^efeh^^), A*^ = -IGDiii^ejelj^i'^). (II.4) 

From now on, we will skip the supperscript "3" for A and A*. Using an appropriate mean field 
approximation, the effective potential of this model can be determined as a function of the condensates 
(a), {tt),{A) and (A*). For simplicity we set (vr) = 0. It is the purpose of this paper to study the effect 



'^ The gauge kinetic term will be added to this Lagrangian in the last step. 
® In Sec. IV and V, the bare mass, mo, will be chosen to be zero. 



of a constant background U{1) magnetic field on the formation of these condensates. To do this, we 
have, principahy, to replace A^ by a classical Aff^ and a dynamical part a^ and then integrate out the 
dynamical gauge field a^ and G^. However, it turns out that for non- vanishing (A, A*), both gauge 
fields A^ and G^ are massive and underlie the Meissner effect/ They are therefore unappropriate to 
be taken as external fields. But, as it is shown in [la . Il9l |. there is indeed a linear combination of 
A^ and G^, that leads to a massless "rotated" C/(l) field, A^ = A^cosO — G^sinO, and a massive 
"rotated" SU{3) field, G^ = A^j^smO + Gf^cosO. According to [l9|, the angle 9 can be determined 
from 

cos9 = :, and sinO = :. (II-5) 

To rotate the fields, one uses the identity 

eQA^ + gT^G^^ = eQA^+~gfG% (II.6) 

and insert the combination OO = 1 on the right hand side (r.h.s.) of this identity. Here, O is an 
appropriate rotation matrix including sine and cosine of 9. The identity pi.6P not only determines the 
new rotated fields as a linear combination of the original non-rotated ones, it also fixes the relation 
between the rotated and non-rotated couplings as eQ = eQ cos 9 — gT^ sin 9, as well as gT = eQ sin 9 + 
gT^ cos9. In the rotated system, one chooses Q so that Q{A) = 0. This leads to 

Q = Q/®lc-l/«>(^)c. (n.7) 

v3 

The above relations between the rotated and non-rotated generators, {Q,T^) and {Q,T), lead then to 
T(A) 7^ 0, which then yields a non-vanishing mass for G^. Hence, as long as the diquark condensate 
A is non-vanishing, the rotated G^ is massive because of T'(A) ^ 0. In this case, the rotated system is 
the true physical system. Once A = and o" 7^ 0, the rotated and non-rotated systems are equivalent, 
because the identity Q{(7) = T{a) = holds automatically [see footnote 8]. Using pi.7p and the above 
relation eQ = eQcos9 — gT^ sm9 between the rotated eQ and the non-rotated eQ, it turns out that 
e = ecos^, as in the electroweak Standard Model. ^ In the six dimensional fiavor-color representation, 
{ur,Ug,Uh, dr, dg,di)), the rotated Q charges of different quarks, in units of e, are presented in Table I. 



quarks 



Ub dr dg db 



A +1 1 _i _i 
2 ^2 22" 



TABLE I: Q charges of quarks in 2SC model in the presence of rotated magnetic field B in units of e. 

Plugging pi.6P in pi.2p . the resulting transformed Lagrangian density is then given by ()II.2p with 
eQA^ + gT^G^ replaced by eQA^ + gTG^, and tt = [see 01. 6h ]. and reads 



Cf = ^Pix)[^J''{^^-^eQA^-igfGl)+fl-f'']i>ix)-mi>{x)^P{x) 

-^A*{ii]^{x)efeh5Hx)) - ^A{up{x)efeh5iP^{x)) - ^ - i|L 



'' As it turns out a is invariant under Uv{i) and SUv{S) groups. Thus Q{a) — T'^{a) = 0, whereas Q(A) 7^ as well as 

T«(A) ^ 0. 
* In a system including mesons and diquarks, only diquarks play the role of electroweak Higgs field. 



with |Ap = AA*. To introduce the external rotated U{1) magnetic field in the third direction, we 
replace A^ — )• yl^^* + a^, with the external rotated electromagnetic field ^^^* in the Landau gauge 
A'^'' = {0,0,Bx,0), and integrate out the remaining dynamical rotated fields a^j and G^. We arrive 
therefore at the full modified bosonized Lagrangian C = Ck + Cf, with^ 

and 

£j = '^{x)[i-f''{d^ - ieQA^^^^) + pL'y^]'ip{x) - Tmp{x)'tp{x) 

-^A*{ilp^{x)efeh5Hx)) - ^A(#(x)e/eh5V'''(^)), (H.IO) 

in a constant (rotated) background U(l) magnetic field B = -663. In what follows, we will simplify 
pi.iop using the method presented in 20] and arrive at an equivalent Lagrangian, which will then be 
used in Sec. Ill to determine the effective potential of the above model in the presence of a rotated 
background U{1) magnetic field B. To do this, we introduce the rotated-charge projectors Jl^, that 
satisfy the eigenvalue equation Qi^g = qVlq. They are given by 

n^i = diag(l, 1,0, 0,0,0) 
^_i = diag(0, 0,0, 1,1,0), 





Oo = 


= diag(0, 0,0, 0,0,1), 




Oi = 


= diag(0, 0,1, 0,0,0), 


and satisfy 







y rjq = 1, and QqQqi = Sqqi. {11.11) 

ge{o,i,±|} 

Using the definition ipq{x) = Qqtl^^x), the fermion field in the six dimensional color-flavor representation 
can now given by 

i;= Yl ^9- (11.12) 

ge{0,l,±i} 

Introducing, at this stage, the Nambu-Gorkov bispinor wave function 

the part of the Lagrangian which is bilinear in ■0, i.e. Cf from pi.lOp . can be brought in the following 
form: 

^f = l E ^5-(^)'5,-*9-(^), (11.13) 

<?G{0,l,±i} 



Comparing to (|II.8|) . in (|II.9|) . we have added the kinetic term of the rotated (7(1) external gauge field ~j{F^,^) 



SI 

2 



where Sq for q £ {0, 1} is given by 



■Sqe{o,i} 



and for q G { — ^, +^} by 



^m-h^+i} 



[S)]-^ 






(11.14) 



(11.15) 



Jj,ab ^ ,-A^«i\afe^,5 ,, „,^i, ,, ^/ = ^,„^t. 

iA*r2A27^. They can be read from pi.iop and the relations pi.lip as well as the definition of Tpq = Q.qil). 



Here, [G% ^-' 



(q)j = b^iid^ + eqA^ ± fi6^o) - m], and k^'^° = iAr^^Xfj^i^ as well as K = 7ok'7o 



Note that pi.lSp can be equivalently expressed as 

Cf = I Yl {??(^)[G(l,-)]-Vf(x)+?,K^)[Gj-)]-V,-(x)-?f(x)4^,-(:E)-V^^-(x)K,-V'f(x)|, 
ge{o,i,±i} 

(11.16) 

where, kg is defined by kg = Qq'K^q. In pi.l6p . Kg is non-vanishing only for q' + q = with q' ^ q. 
|, +^} we have therefore 

' +iA^^ if {p,a) = (2,4), 

jp^- ^ -iAj^ if {p,a) = (1,5), (11.17) 

otherwise, 



For q £ {"li+^l we have therefore 



(«^5=_i)p<7 = (f^i«:f^_i)pff = (/^f^.i) 



l'^q=+iJp<T 



(ii_iKrii )„o- 

2 2 



(kO_ 



ypcr 



+iAj^ if (p,a) = (4,2), 

-iAj^ if (p,a) = (5,1), (11.18) 

otherwise, 



whereas for q G {0, 1}, we have 



Kij = (S7jK%)p^ = 0, for i,j = 0,1 or i = j = ±-. 



(11.19) 



This is in contrast to the case of three-flavor color-flavor locked (CFL) phase, studied in 20|. In that 
case, there exists a charge q = —1 and the combination of (ri=piKri-|-i) leads also to nonzero result. 



III. ONE-LOOP EFFECTIVE ACTION AND THERMODYNAMIC POTENTIAL 



In what follows, the one-loop effective action of the theory, F, will be determined in the mean field 
approximation in terms of o" = {a{x)),A = (A(x)), and A* = (A*(x)). Using the following path 
integral over the quark fields 



ir[o-,A-A* 



Vil^Vil; exp [ i d x C 



(III.l) 



where, C = Ck + Cf, with C^ and Cf from pi.9P and pi.l3p . the effective action up to one- loop 
quantum corrections is given by 

.2 1^12 q2^ 



r[a,A,A*] 



O" 



AGs AGd 2 



+ i^ + ^)i^ + riff^[^'^'^ 



(III.2) 



Here, V is the 4-dimensional space-time volume, and r),^ is the one-loop contribution to the effective 
potential. It arises by integrating out the fermion fields and reads 

riff = -^ E Tr{7VGc/sx} HSf% (ni.3) 

where Sq is defined in pi.l4p and (jll.isp . Here, the trace "TV" operation in (jHLSP includes apart 
from a two-dimensional trace in the Nambu-Gorkov (NG) space, a trace over the whole phase space. 
It is therefore defined by a trace over the color (c), fiavor ff), and spinor (s) degrees of freedom, as 
well as over a four-dimensional space-time coordinate (x) [g]. To compute pil.sp . we have to notice 
that, according to Table HI the blue quarks {ub,dh) have q = 0,1, whereas the green and red quarks 
{ur , Ug , dr , dg) have q = it^. Thus relation (|III.3p reduces to 

ri^^= E ^«/^ (III.4) 

where the one-loop effective action of the blue (b) and red/green (r/g) are given by 



^ riff^/" = --Tr|;vGc/sa}ln[cS_^i"^] - -Tr|jvGc/sx}ln[5_i"^]. (III.5) 

Kg{r,g} 

To perform the trace operation in the NG space, we use 

det ( "^ ^ j = det {-BC + BDB-^A) = det {-CB + CAC-^D) . (in.6) 

Using further tr In j4 = In det A, we arrive at 

detjATGc/sxli^o""^] = detjscx} [{7"(«5a + fiSao) - m}{7"(i5a - ix6ao) - m}] , 
'^^^{NGcfsx}[S+l^^\ = detjscx} [{I'^iida + eAa + fj-Sao) - 'm}{j°'{ida + eAa - jlSao) - m}] , 
det{jvGc/sx}['5±i"^] = det{,^^y[\A\'^ + {--f'^{ido,±-eAa + fl6ao)-m} 

x{7"(ia, ± hAa - fi6ao) - m}] , {111.7) 

where we have skipped the superscript "ext" on the external rotated gauge field A^. Here, fl = /x — 2/^8; 
fL = fj, + fis and m = niQ + a. The determinants in pil.7p are now to be calculated in the momentum 
space. To do this, a generalization of the method described in [2Q] for arbitrary charges is necessary. 



This method is originally developed by Ritus in [27| in order to determine the Green's function of 



charged fermions in the presence of background magnetic field. It is then extended to charged vector 
fields in [2g] . Recently, it is used in [29] to determine the electric-current susceptibility of quark matter 



in the presence of external constant magnetic field. As it is described in 29|], in the Landau gauge for 
the external rotated gauge field, a projection operator P„ can be defined 



Pn = kifn^x) + /„_(X)] + i[/n+(x) - fn.i^Wl^ ^OV qB > 0, 

Pn = UfnM + fn-{x)] " i[/n+(x) " LA^Wl^ fo^ qB < 0, 



(III.^ 



that includes the basis functions /„^ (x) defined by 



fn+{x) = 0n(x-||), n = 0, 1,2,--- , 

fn.{x) = (/)„_i(2;- 1^), n = 1,2,3- •• . 



qB 

Here, (pnix) are the standard Landau quantized wave functions [29] 



(t>n{x) 



-f\qB\Y/' 



2"n! V IT 



] exp (-^\qB\xAHn [V\qB\x) , (HI- 10) 



with Hn{x) the Hermite polynomial of degree n. Using the projectors Pn from (jIII.Sp . it is easy to 
show that 

7/, {idf" + qA^') p^e-'^vot-Pyy-P^z) ^ p^ fo _ sgn(qB)^2\qB\n-f^ - p,-f^^ (.-i{pot-Pyy-P.z) ^ 

(III.ll) 

The r.h.s. of pil.lip is a free Dirac operator with a modified momentum 

p^ = (po,0,sgn(5S)V2kS|n,p3)- (111.12) 

This shows also that the solution of the Dirac equation in the presence of a constant magnetic field 
can be given by a combination of the projection operators Pn and the ordinary free Dirac spinors 
u{p,s) and v(j),s) (291]. 

To compute the determinants in (jIII.Tp in the momentum space, we will use, for the charges q ^ 0, 
an appropriate momentum basis, similar to piI.12p . and for q = 0, the ordinary four-momentum p^. 
In other words, we have 

P'^^o = iP0,0,jli^2\qeB\n,p3), for q = l,±i, (m 13) 

Pf=0 = iPo,Pi,P2,P3), for q = 0, 

where A replaces sgn(qB) in piI.12p . This leads to the well-known quasiparticle dispersion relations 



k1 
in the presence of a constant magnetic field aligned in the third direction [23l ]. 

Eg = ^2\qeB\n + pi+w? for q = l,±i, 

£^0 = -s/pi +P2+P3+ m"^ for q = 0. 

Using the momenta (JTTLTSJ) and transforming pil.4p and pil.5p into the Fourier space, the one-loop 
effective action reads 

fi«^(p)= E fi'J^'^iP), (111.15) 

KG{r,g,fe} 



with 



fi«^^'(p) = -i Y. i^detMEq+^^f-pl}{{Eg-flf-pl}], 
ge{o,i} 

Y, ri'J^'-ip) = -2i Yl lndet..[(4+)'-p2)(i?(-)'_p2)]. („i.i6) 



Here, Eg for q G {0, 1, ±i} are defined in (|III.14p . and E^^'^ = ^ {Eg ± fif + |A|^ for q G {+i, -i}. 
The factor 2 in tlie last equation of pil.iop reflects the degeneracy in the quark charges for Uj./g as 
well as dr/g (see Table H]). Note that a trace over Landau levels, n, is implemented in the expression 
on the r.h.s. of piI.16p . This is because E^^q from piI.14p depends explicitly on n. This trace will 
be performed in the next step, where the one-loop effective action will be explicitly determined in 
the momentum space. Performing the remaining determinant in the coordinate space leads, for a 
constant background magnetic field, to a space-time volume V. At this stage, we will introduce the 
effective thermodynamic (mean field) potential Q^^, that is defined by the effective action through the 
relation Fefi = — Vfioff- To determine the one-loop contribution to the one-loop effective potential at 
zero temperature 0,],J, it is convenient to determine it first at finite temperature, and then taking the 



limit T — )• 0, consider only the zero temperature effects 23|]. For quarks with q = 0, one replaces 
Po by po = icoi,^^ where uj£ are the Matsubara frequency defined by w^ = {2i + l)vrT, and the po 
integration by an infinite sum over the Matsubara frequencies. For an arbitrary function f{po, pq), we 
get therefore 



f d'^p 1 ^ /■ d^p 

J (^^(^0'P5'=o) = ^ z2 J (^ fi^^i^P)^ (111.17) 



where P = T ^ is the inverse of the temperature T. For the quarks with q ^Q, apart from a summation 
over the Matsubara frequencies (., a summation over the Landau levels n is also to be considered [see 



(1111.131) ]. We get therefore [36 1 

/ 7^:^/(Po,P9Vo) = ^— 2^ 2^an / -^ f[iuJe,n,P3), (IILIS) 



where a„ = 2 — 5no refiects the fact that Landau levels with n > are doubly degenerate 20|, [23 1 . 



Following the above recipe, the one-loop contribution to the thermodynamic potential is given by 

9e{o,i,±|} 
where for q = 0, we have 

J^f=o = iV-^lndet[(3^{{Eo + ilf + ujj}{{Eo-ilf+ioj}] 

= -^ E J^H^'{iEo + iif+u:j}{iEo-flf + u:j}] 

= --pj -§^, [PEo + In (l + e-/^(^o+^)) + In (l + e-/^(^o-^)) } , (IIL20) 



^^ The effect of tlie cliemical potential is already considered in rj,jj as well as Qr^g . 
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and for q = 1, we have 



+00 +00 +00 J 

, n J— 00 "^ 



^= — 00 71=0 



+00 f+OO J 

Y^anj g {pE+1 + In (1 + e-^(^+^+«) + In (l + e-''(^+-^)) } (III.21) 



/3 
Note that £"+1 depends exphcitly on n that labels the Landau levels. Finally, for q = iti, we arrive at 



Y, ^L = 4iV-i Indet [/3'^{S+ ' + coJHE.^ + c^|}] 



9e{+i-^} 



hOO +00 /■ + 00 



, n J — 00 J-DTT 2 2 



^ e=-oo n=0 

+°o /.+00 



(111.22) 

where S^i = £ti and E^-, = EZ^ are used. Plugging pil.20p - pil.22p in piI.19p and taking the limit 



22 22 

r — )• by making use of the relation [37 1 



lim Tin [1 + e'r] = -xOi-x), (III.23) 

with 9{x) is the Heaviside 6'-function, the temperature independent part of the effective potential, 
including the tree level and the one-loop corrections reads 

-2 y ^ [i?o - {Eo - /i)^(/i - Eo) - (/i + £;o)^ (-/i - Eo) ] 

-eB y^an /I [£+1 + Et^ + E^, + {fl - E+,)9{fl - E+i) - 2E+ 9{-E+) 

n J-OO 47r^ — — +2 +2 



n=0 



-2i?7i0(-i^7O-(/i + i?+i)0(-/i-i5;+i)]. (111.24) 



The above result pil.24p is comparable with the result in [6|], which is derived for a similar 2SC model 
in the absence of the magnetic field B. In this case the thermodynamic potential up to one- loop order 
at finite T is given by 

K£{r,g,fe} 

where for different colors, we have 

f^W/* = -2iV"ilndet[(^2 +^2^(^2 ^^2)] 

= -^l^{PE + In (1 + e-^(^+^)) + In (l + e-'^(^-^)) } , (III.26) 
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and 

Ke{r,g} 



y"^{/3i?^ + /?i?^ + 21n(l + e-^^-)+21n(l + e-^^t)},(III.27) 



£=-oo 

_ 4 /■ d^p 

with E = v^p2 + m2, E± = E ± fl, and S^ ee y^(S ± /x)^ + |A|2. Using ()III.23p . the temperature 
independent part of pil.25p reads 

n-2 IA|2 r rfin 

-(/i + E)0(-/i-E)]. (III.28) 

IV. ANALYTICAL SOLUTIONS OF THE xSB AND CSC GAP EQUATIONS IN THE LLL 
APPROXIMATION: A COMPARISON OF B = AND B^Q CASES 

In the previous section, the one-loop effective action of the NJL model including meson [a) and 
diquark (A) condensates in the 2SC phase at finite ei?, /U and ^s is computed in the mean field 
approximation. This is the purpose of this paper to have a complete understanding on the effect of 
external magnetic field on the formation of these condensates. To this purpose one has to solve the 
following gap equations and color neutrality conditions 



da 



0, 



dA 
dQ,B{a,A,H8;iJ',eB) 



8^8 



o-bAb,ij.s 



0. (IV.l) 



The solutions of the first two equations build the "local" minima of the theory. In Sec. IV, we 
will solve the above equations numerically for any value of the rotated magnetic field eB. Keeping 
{a^ A) 7^ (0, 0) and looking for global minima for the system described by complete JlcffCo", A, /Ug; /x, eB) 
from pil.24p in the presence of the rotated field, it turns out that in the regime 300 ^ /U < 500 MeV, 
the system exhibits two "global" minima. They are given by {(Tb t^ 0, A^ = 0, /ig = 0) in the regime 
/i < ^c-: and [as = 0, A^ 7^ 0,/i8 7^ 0) in the regime /i > /ic- Here, //c is a certain critical chemical 
potential, and, shall be determined numerically in Sec. IV for a wide range of eB [see Fig. 9]. We will 
denote the regime characterized by {as / 0, A^ = 0,/i8 = 0) and {ub = 0, A^ / 0,;U8 / 0), by the 
xSB and the CSC phases, respectively. In this section, we will analytically determine the solutions of 
the above gap equations in the limit of strong magnetic fields \qeB\ ^ /x^, and in the xSB and the 
CSC phases separately. We will then compare these solutions with the corresponding solutions of the 
gap equations in i3 = case. In the above limit, the dynamics of the system is dominated by LLL. 
The goal is to determine analytically the mass gaps of the xSB and CSC phases separately. This will 
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be done in Sec. IV. A and IV. B, respectively. In IV. A. 1 as well as IV.B.l, we consider the case of 
strong magnetic field, whereas IV. A. 2 as well as IV. B. 2 are devoted to i? = case. 

A. The chiral symmetry breaking phase 

1. Strong magnetic field 

According to the descriptions from the previous paragraph, the xSB phase is characterized by ctb t^ 
and As = /ig = 0. To study this phase in the LLL approximation, we will, in particular, focus on the 
first gap equation from pV.ip 



or equivalently on 



11 



aa'r(^,A,/i8;^,e5) 



da 



0, 



(IV.2) 



CTs,Afl=/i8=0 



aOLLL(^,AB = ^8 = 0;^,e5) 



da 



(IV.3) 



where J7^sb arises from pil.24p with n = and {as / 0, A^ = ^g = 0). To solve pV.SP analytically, 
let us consider ri^gg first in the momentum space 



n^lia, Ab = /X8 = 0) = ^ + ^ - 2 



(2^ 



\/p2 + 0-2 + [^- Y^p2 _^ o-2)6'(^ - A/p2 + 0-2) 



2,eB 



a' B_ 

4Gs^ 2 



-3eS 



47r2 



pl + a^ + {tJi- Jpl + a^)e{^ - Jpl + a 



\ '-^^/p^ + a^ + e{^i-a) / 



p^dp 



vr^ 



(/i - Vp^ + a2) 



^V^^T^ + 0{fi-a)l^ ' 



(/"- VP3^ + <7^J 



27r2 



(IV.4) 



Here, we have introduced the momentum cutoff A for the first integral arising from the contribution 
of zero charged particle. In contrast, the momentum cutoff A^ = \/eB is chosen for the first inte- 
gral proportional to eB, that arises from the contribution of the remaining three charged particles. ^^ 
Considering furthermore the effect of the Heaviside ^-functions in the integrations limits, the corre- 
sponding momentum cutoff to the remaining two integrals is given by y^fj,"^ — a^ with the assumption 
that a < fi (see the 9{iJ, — a) before these two integrals). Performing the integrations over p = |p| and 



^^ In [Tj, the same procedure is performed to study the x^B and the CSC phases separately. 

^^ The charges of the particles is defined with respect of the rotated magnetic field. They are presented in Table 1, in 
units of e. 
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P3 in pV-4p . we arrive at 



fi^^^(a,AB = ^8 = 0) 



a^ B' g^ 
A 



In 



A + VA2 + Cj2 



0(/i-CT)ln 



/i+ V/U2 - 0-2 



O" 



g^, 7a^T^(2A2 + ^2) - ^^(;u - a) v^ 



^2 (2/i2 - 5(7^ 



3eScr^ 
2,eB' 



In 



Ar + JA2„+CT 



B 



0(^ - a) In 



;U+ \//u2 - 0-2 



O" 



47r2 



AbJ^I + 0-2 + 6'(^ - (j)n^J^i'^ - 0-2 



Minimizing the above potential according to pV.3p . the gap equation reads 






9a 



+o-|jln 



■K 

G's 



(^B { T^ Ao-ba/A2 + a%- 3eBln 



'B 



Kb + a/A| + al 



o-fi 



A+ JA2 + 



a 



b 



CTB 



+ e{n-aB) 



lx\Ufi -a% + {3eB - a%) In 



fJ-+ \ fJ- 



a 



B 



CTB 



(IV.5) 



> . 



(IV.6) 



To find a nontrivial solution ctb 7^ to this equation, we expand it in the orders of the dimensionless 
and small parameter x = ^ <C 1 up to order O(x^), and get 



2AM - 



2 /'S , /4A 



B 



+ 3eB In ' ^ 



o-f 



+6'(/i - aB) { 0% - 2^' + {a% - 3eB) In ( %- 



a 



B 



(IV.7) 



where the dimensionless coupling gs = ^-2 is introduced. In what follows, we consider two different 
regimes of // < as and ^> gb separately. To find real solution for the simplified gap equation pV.7p , 
we will then distinguish various regions for the dimensionless coupling gg. 
i) In the first regime characterized by /i < crs, the gap equation pV.7p reads 



2A2 1 



5^2 2 1 /^4A2\ _^^ /4A2„ 



o-filn 



o-r 



(IV.. 



Since for < 5^ < 1, the l.h.s. of PV.SP is positive, a nontrivial real solution arises only by the 
assumption a^ln (^4-) ^ 3ei?ln ( ^5^ ) , which is indeed justified in the LLL approximation. Ne- 
glecting therefore the first two terms on the r.h.s. of pV.Sp . we arrive at 

o-| = AeBe ^-^^Vs. V. (IV.9) 

Note that the assumption a^ In f ^4- ) ^ 3e5 In ( —3^ j does not set any limitation on the relation 
between two momentum cutoffs A and A^. Depending on whether A is larger or smaller than A^, 
different regimes are to be distinguished for the coupling gs'. 

(IV. 10) 



For A < Kb-, we get < 5^ < — 



A2 



For A > As, we get < 5^ < ^^^gig^^^ < ^ 
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The dynamical mass as from pV.9p is, apart from numerical factors, the same as the dynamical mass 



of the NJL model in the presence of constant magnetic field from 3l[. The additional factor 1/3, that 
arises in the exponent of pV.9P corresponds to three different quark charges q = 1,±2 that have, in 
the regime of LLL dominance (n = 0) equal contributions to the effective potential in the xSB phase. 
Let us consider again the gap equation (|IV.8|1 for the case Qs > 1. In this case a nontrivial solution 
may exist only for A in the same order of magnitude as A^. To find the solution, we rewrite first the 
gap equation (jIV.SP as 



\2 



Expanding now the second term on the r.h.s. in the orders of y = ^^ — 1 ~ up to 0{y'^), we arrive 
at 



3K 



2AM^-l|»-,|ln^g-3a(||-l)+3.B.n(^|. (IV.12) 



Using the same approximation a^ In ( ^3- I <C 3eB In I — j^ I , we can neglect the first term on the r.h.s. 
of (|IV.12p . and arrive at 



+ 3eB, (IV.13) 

whose solution reads 

a% = CA^e'^^isii, with C = -e ^ 7.25. (IV.14) 

ii) In the regime characterized by as < A*, the gap equation is given by 

It arises by expanding pV.6p in the orders of x = ^ up to order 0{x^). As it turns out a real solution 
may be found by expanding pV.lSp in the orders it; = ^ — 1 ~ and z = -^ — 1^0 up to 0{'uP') 



B 



as well as 0{z ). The mass gap can be computed directly from the resulting equation and reads 

2 ,2/1 \ 3e5 5u2 ^^^^ ^ 

Note that a real solution for gb in this regimes arises only when eB and Qs satisfy the following 
conditions: 

/ M^ < Ab < ^, and 3g^4A'_3 o < 5s < 1, 
For < 5, < 1 : <^ 5 2 ^ ^ 2!^ ^ oA^ 

[-1<^B, ^ and 3,-b41._3^2 <9^s< ,,B+li-._,^. , (IV.17) 

For g, >1: \^^^<A^<^^ and K g. < 3,-^41' ,5^. • 

The above conditions arise without specifying any relation between A and A^. 
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2. Zero magnetic field 



To clarify the effect of strong magnetic fields on the mass gap, we will present in what follows the 
analytical results of the gap equation corresponding to the effective potential pil.28p at zero magnetic 
field. ^^ Setting, as in the previous section, in the xSB phase, Aq = /^g = in the corresponding 
effective potential ()III.28p . the resulting potential in the momentum space reads 

.2 M p2^p 



5^xSb(o-,Ao = ^8 = 0) 



AGs 



6 



^" 



(IV. 18) 



After performing the p integration, we arrive at 

-2 3m^ 



^xsb(o",Ao = /Ug = 0) 



3A 



a 



+ 



4Gs 47r2 



,„fA±^^5!±Z _,(„_, ,iJfi±VZZZ 



a 



a 



-^(2A2 + a^)^A^ + a^ + -f^{5a^ - 2/x2)e(/i - a)v^^I^^^. 



(IV. 19) 



The corresponding gap equation reads then 







'xSB 



da 



^-AyX^.,.^,„ A±4^ 



o-o 



+e{fi - ao) 



fiJu^-a^-aQln 



/i + ^//? 



CTn 



o-Q 



(IV.20) 



Defining, similar to what we did in the B ^ case, a dimensionless parameter x = ^ ^ 1, and 
expanding the gap equation pv.20p in the orders of x up to 0{x'^), we arrive at 



an 



XSB 



da 



2A^ 



o-Q ( 1 - In —2- ] +9{fi- CTo) ( -o-Q + 2//^ - CTq In -^- 



(IV.21) 

where we have introduced the dimensionless coupling Qs = — t — = Qqs- To find a real solution for 
the mass gap ctq, we have to distinguish two different regimes of the chemical potential /i. They are 
characterized by ;U < o"o and do < /i. 
i) In the first regime characterized by /U < (Tq, a real solution (Tq < A arises only for 1 < ^^ < 5i with 



^1 = Y^T^ — 3.26. It reads 



al =4A2exp( W^i 



1 

Ws 



(IV.22) 



It corresponds to one of the two real branches of the Lambert M^(x) function, Wq{x) and W-i{x), 

38 1 for more details on the Lambert VF-function)^^. 



which is known to be the function satisfying (see 

W{x)e^^''^ 



X. 



(IV.23) 



^^ In [6(1 the full gap equations of the 2SC model including the mesons is solved numerically for B = 0. As it turns out 
the system exhibits, as in i? 7^ case, a phase transition from the xSB to the CSC phase. Here, the xSB phase is 
characterized by (ctq 7^ 0, Aq = 0) and the CSC by (cro = 0, Aq 7^ 0). 

^* According to the explanations in [3g]: If x is real, then for — 1/e < a: < 0, there are two possible real values of Wix). 
One denotes the branch satisfying — 1 < W{x) by Wo (a:), and the branch satisfying W(x) < — 1 by W-i{x). 
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a) In the second regime characterized by (Tq < /i, the gap equation is 



AM - - 1 
.9s 



2 2 , ^0 1 / ^ 

0"n - ^ + ^ In ( — 



(IV. 24) 



It arises from pv.20p by an expansion in the orders of x = ^ up to 0{x^). Introducing a smaU 
parameter w = -j^ — 1 and expanding pV.24p in the orders of w up to 0{'w'^), yields the mass gap 



a^n = 2A2 



9s 



1 +2fi\ 



Note that a real solution for (Tq arises only for ^^ > 1, satisfying 



2A2 



A2 



-2<9s< 



,2' 



(IV.25) 



(IV.26) 



2A2 - /i2 ^ ^^ ^ A2 - /x- 

In Sec. V, we will perform a numerical analysis to study the xSB phase for any arbitrary magnetic 
field. We will show that, similar to the B ^ case, the second regime characterized by do < /i belongs 
to the color symmetry breaking phase and is indeed irrelevant for the present xSB phase. Comparing 
therefore only the relevant part of the solutions, i.e. pV.9p for B ^ with ()IV.22p for i? = 0, we 
note that, in contrast to B = case, in the presence of strong magnetic fields, the formation of chiral 
symmetry breaking bound state as is possible even for small dimensionless xSB coupling < ^^ < 1- 
This is in fact one of the consequences of the phenomenon of magnetic catalysis 3l| . 



B. The color superconducting phase 



1. Strong magnetic field 



According to our explanation in the paragraph below OV.lh . the CSC phase is characterized by 
(as = 0, A 7^ 0, /ig 7^ 0). The corresponding effective potential arises from pil.24p by setting n = 0, 
and as = 0. In the momentum space, it is given by 

'0 

'■°°dp3 



^csc(^i?=0,A,/X8) 



A2 ^2 ^oo 



b + (A - pWifj- -p)-{fj- + pW {-P' - p)] 



-eB 



27r2 



Ps + iP-- PsWifj- - Ps) - (/i + P3)0 (-/i - Ps) 



+ViP3 + fiy + A2 + v(P3 - fiy + A2 



(IV.27) 



Performing the ps and p integrations by introducing the momentum cutoffs A and A^ for the p = |p| 
as well as pa integrations,^^ we arrive at 



n^fci^B = 0,A,fis) 



127r2 



A2 B^ 3A4 + {3eB + il'^)fl^ 
4G^^T 
eB 

~4^ 



+A^ln 



A| + {Ab - /I)VA2 + (Ab-m)' + (Ab + /2)VA2 + (Ab + /2)2 

(IV.28) 



2 ^ , VA2 + {Ab + /i)2 + (Ab + ft) 



VA2 + {Ab - fiY - {Ab - n) 



See our description below (|IV.4[I for the choice of the momentum cutoffs. 
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In the CSC phase, we have a set of two coupled equations: the color neutrality condition, S^'^ = 



anpLL 



9^8 

and the gap equation "'g^'^ = 0. The goal is to solve these two equations to determine fig and A as a 
function of the external magnetic field B and the chemical potential /i. Let us first consider the color 
neutrality condition 



_ OSZ^gQ 



dfis 



Ajj-' + 3eB{2jl+JAj, + {Ab - fiy - J Aj, + (A^ + /i)^ , (IV.29) 



where p, = fi + fig and fl = ^ — 2 fig. Defining three dimensionless (small) parameters, x = ^, y = —^ 
as well as z = ^, and expanding ()IV.29p in the orders of x,y and z up to 0{x'^),0{y'^) and 0{z^), 
/ig can be determined from the resulting equation and is given by 

As for the gap equation corresponding to A^, we have 



n -"CSC 



dA 



7r2 / jAl + {K+-^f + {KB + p) 
— - eB In ) I . (IV.31) 



After expanding ()IV.3ip in the orders of y = -^ and z = ^ up to 0{y^) and 0{z'^), and replacing 
/ug from ()IV.30P in the resulting equation, we arrive at 



+ 7^-eBln{^—^^^]=0, (IV.32) 

gd 9eB V A| y ' 

where the dimensionless diquark coupling in the CSC phase gd = — S— is introduced. The diquark 
mass gap A^ can then be determined directly from pv.32p and reads 

Al = 4(A| - f,^) exp (-^^) ■ (IV.33) 

This result is comparable with the results by [20] for the three- flavor CFL model. In particular, in 
both models the exponents are proportional to {eBg^)^^. The dependence of A^ on the magnetic 
field demonstrates the effect of magnetic catalysis J3l|, that states that even for small value of the 
dimensionless diquark coupling g^, the presence of a strong magnetic field leads to color symmetry 
breaking and the formation of diquark mass A^. 



2. Zero magnetic field 

We consider, as next, the effective potential of the 2SC model in the absence of magnetic field from 
pil.28p in the CSC phase by setting (cro = 0, Aq 7^ 0, fig ^ 0). In the momentum space, the resulting 
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potential is then given by 



^csc(o"o = 0, A,/i8) 



A^ 



A J2 



p dp 



^Gd Jo tt 

+(/i - p)e{jl -p)-{fl+ p)9 (-/i - p) 



VA2 + {p + /2)2 + y/A2 + (p-p,y+p 



A2 /i4 ^4 1 



i^ - TTT - TTrT[3A(A' + 2A2) + ^(ISA^ - 2A2) - 2/i2(A + /2)] y^A^ + (A - A)^ 



AGd Qir"^ 27r2 127r2 

+ [3A(A2 + 2A2) - /i(13A2 - 2A2) - 2fl^{A - /2)] ^A^ + (A + /i)2 



3A^(A^-4/z2)ln 



■ [A + /2 + VA2 + (A + /Z)2] [A - /2 + VA2 + (A - /i)2] ^ 



A2 



(IV.34) 



In this case, the color neutrality condition reads 

C*' 'CSC 







5/i8 



3 .A2.1 / [A + /i+VA^ + (A + /2)2][A-/2+VA^ + (A-/2)2] 



2/i-*-3A^/xln 



Ao.As 



A2 



[2A^ - A^ - fiiA + fi)]JAl + (A - /2)2 + [2A^ - A^ + /i(A - /i)] JA2 + (A + /Z) 



(IV.35) 



After expanding in the orders of x = ^, y 
0{z^), we arrive at 



/^8 



-^, and z = ^ up to order O {x^), 0{y^), as well as 

(IV.36) 



A 2 A 2 / A 2 



3// 6/i V4A2y' 



where Aq satisfies the gap equation 



a 



CSC 



dA 



AcAs 



IT 

'2G^ 



+ (A + 3/z) VA2 + (A - /i)2 + (A - 3/Z) Ja2 + (A + /x)2 



,(^2 _ 2^2)1^ / [A + A + VA2 + (A + /x)2][A-/2+VA2 + (A-A)2] 



(IV.37) 



Using the same method as above and expanding pV.37p in the orders of y and z up to order 0{y^) 
as well as 0{z^), we get 







Q.r 



dA 



1 



Ao,A8 V 5dy V4(A2-/i2); 



A2 



(IV.38) 



where gd = — % — = Ag^- Solving pV.38p . the diquark mass for vanishing magnetic fields is then given 

by 



Al = Ca (A2 - ^2) exp (-^(^-l\\ with Cs = 4e'3 ^ q_2. 



(IV.39) 



The qualitative behavior of Aq as a function of fi coincides with the results from [6|, |39[. The color 
chemical potential jls arises by replacing (1IV.39P in (IIV.36J) . These results are to be compared with 
PV.33P [for the diquark mass gap A^] as well as ()IV.30P [for the color chemical potential fis] in the 
presence of strong magnetic field. 
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V. NUMERICAL RESULTS FOR ARBITRARY MAGNETIC FIELD 

In the previous section, we have presented analytical solutions for the order parameters a and A 
corresponding to x^B and CSC phases, as well as for the color chemical potential /ig in the presence 
of strong magnetic fields in the LLL approximation. We have then compared our results with the mass 
gaps arising from the thermodynamic potential of the 2SC model in the absence of magnetic field in 
order to emphasize the effect of strong magnetic fields on the formation of bound states a and A in 
the superconducting 2SC model. In this section, we will study numerically the effect of any arbitrary 
magnetic field on quark matter without restricting ourselves to LLL approximation. In particular, we 
are interested on the dependence of the mass gaps on the external magnetic field eB and the chemical 
potential /x. To do this, we set, as in the previous section, niQ = and choose Gd < Gs- Comparing 
our numerical results with the analytical results from Sec. IV, we will determine numerically the range 
of the magnetic field strength for which the LLL approximation is reliable. At the end of this section, 
we will study the phase diagram of the model in a ;Uc — eB plane, and determine the type of various 
phase transition between the xSB and the CSC phases for a wide range of eB. 

Let us start with the one-loop effective potential ()III.24p arising from a mean field approximation 
in the presence of an arbitrary magnetic field. To perform the momentum integrations numerically, 
we have to fix the free parameters of the model, the momentum cutoff A and the couplings Gs and 
Gd- Our specific choice of the parameters is [91] 

A = 0.6533 GeV, Gs = 5.0163 GeV-^, and Gd = jGs- (V.l) 

For vanishing magnetic field eB = 0, they yield the xSB gap (Tq ^ 323.8 MeV at /x = 250 MeV, and 
the 2SC gap of Aq ~ 126 MeV at /i = 460 MeV.^^ Smooth cutoff functions (form factor) 

l + exp(M^) l + exp(^^2±ap^) 

are then introduced to perform numerically the momentum p integrations corresponding to zero 
charged particles and charged particles, respectively.^ In ()V.3p . ^ is a free parameter and is chosen 
to be ^4 = 0.05A. Similar smooth cutoff function (form factor) is also used in |22l ]. Here, as in |22l |. the 
free parameter A determines the sharpness of the cutoff scheme. In what follows, we will first study 
the behavior of mass gaps as well as magnetizations in the xSB and CSC regimes as functions of eB 
and for fixed chemical potentials. 

In Fig. la, the xSB mass gap a is plotted as a function of eB and for fixed chemical potential /i = 250 
MeV. Small oscillations for small value of eB arise from the well-known van Alfven-de Haas (vAdH) 
effect 2j]. They occur when the Landau levels pass the quark Fermi surface. They are also observed 
in 35| for the xSB mass gaps. Note that the oscillations are sharper, the smaller the value of the 



^^ Although our free parameters A, Gd , and Gs coincide with the parameters used in [9| , the numerical value of cro is 

different from what is reported in [9|. The reason for this difference is apparently in the choice of the cutoff function. 

Whereas we use smooth cutoff function (|V.2p . in [9| a sharp momentum cutoff is used to perform the momentum 

integrations numerically. 
^^ In HIII.24[) . the integrals proportional to eB and including a summation over Landau levels n arises from charged 

quarks with charges q — ±^,+1. 
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FIG. 1: a) The dependence of ctb on eB in the xSB phase for fj, = 250 MeV; b) The analytical result of (Tb in 
the regime of LLL dominance from (|IV.9p is plotted for eB G {0.45, 1} GeV^ and fj, = 250 MeV. 



free parameter A in (jV.SP is chosen [see also 40|] for a discussion on the effect of free parameters in 
smooth cutoff functions (form factors)]. ^^ As for eB > 0.45 GeV^, where the dependence of aB on 
eB is almost linear, we enter in the regime of LLL dominance. The qualitative behavior of as as 
a function of eB for strong magnetic fields can be checked by comparing our numerical result from 
Fig. la with the analytical result for as from pv.9p .^^ The latter is plotted in Fig. lb for the same 
interval of the magnetic field, i.e. eB G {0.45, 1} GeV . Similarly, in Fig. 2a, the CSC mass gap 
A^ is plotted as a function of eB for ^ = 460 MeV. Same small vAdH oscillations appear for small 
eB < 0.47 GeV^. They are also observed in 2^] and 23|] for the diquark in the CFL superconducting 
phase. Small oscillations in Fig. 2a, end up in a linear regime, that starts, as in the previous case, at 
eB > 0.47 GeV^. The qualitative behavior of A^ in this regime can be compared with the analytical 
result (|IV.33p . that arises in the LLL approximation (Fig. 2b). 
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^* We have also checked our results for A = O.OOIA (quasi-sharp cutoff), where instead of small oscillations, small 

discontinuities appear in the regime eB < 0.4 GeV^. 
^^ For our specific choice of Gs and A from (|V.H) . the dimensionless coupling < Qs < 1. On the other hand, since no 

mixed phase is assumed here, as from pV.9|l in the regime ^i < as is the only relevant mass gap that can be compared 

with as arising from our numerical calculation. 
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FIG. 2: a) The dependence of A^ on eB in the CSC phase for /i = 460 MeV. b) The analytical result of A^ 
in the regime of LLL dominance from (|IV.33P is plotted for eB e {0.47, 1} GeV . 

In Fig. 3, the dependence of the color chemical potential /ig on eB is plotted for fi = 460 MeV in the 
CSC phase. The vAdH oscillations in Fig. 3 are similar to the oscillations of fis in the regime of small 
magnetic fields that are observed in 2^ in the superconducting CFL model. 
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FIG. 3: The dependence of /is on eB in the CSC phase for /i = 460 MeV. The vAdH oscillations are similar to 
the oscillations that are observed in 



22| in the superconducting CFL model. 



In summary, comparing the above numerical results with our analytical results from Sec. IV. A 
and IV. B for non-zero magnetic field, it turns out that there exists a threshold magnetic field {eB)t ~ 
0.45 — 0.50 GeV^, where the qualitative behavior of our numerical results coincides with the qualitative 
behavior of the analytical results for xSB and CSC mass gaps ub and A^.^'' This regime, for which 
the LLL approximation seems to be reliable, will be denoted from now on by "the linear regime". ^^ 

Using the above data, the magnetization of the 2SC superconducting medium can be studied as a 
function of eB and for fixed chemical potential fi. Fig. 4 shows the dependence of the product of the 



do. 



(1) 



magnetization M = Me^ with M = — q^ and the rotated magnetic field B = Be^^, as a function 
of eB for two different chemical potential ^ = 250 MeV in the xSB regime and // = 460 MeV in the 



)(i) :. 



CSC regime (Figs. 4a and 4b). Here Q^'g' is the one-loop part of the effective potential pil.24p . For 



,9n 



(1) 



simplicity, we will use the definition M-B = —eB qS^ . Equivalently, one can define the magnetization 
by introducing the Gibbs free energy density Q in the presence of a constant magnetic field B 



B' 



,Wi 



g{a,A■,B,^l) = — + nit,' {^, A; S, ^) - HB, 



(V.3) 



where H is the external magnetic field [231 ]. Whereas in vacuum H = B, in a medium with finite 
magnetization density, the external magnetic field H is different from the induced magnetic field 
B. Minimizing Q with respect to B and evaluating the result at the minimum of the potential, we 



^'^ Note that the similarity in tlie numerical and analytical results for eB > {eB)t is only qualitative. This is because of 
various approximations that are carried out to determine the analytical results [see Sec. IV for more details]. 
Note that in Figs. 2 and 3, the threshold magnetic field satisfies the requirement of LLL approximation {eB)t ^ fj- . 
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get the well-known relation M = B — H, where M is the magnetization. Note that the minimum 
of the potential in the xSB phase is given by {as 7^ 0, A^ = ^ug = 0) and in the CSC phase by 
{as = 0, As 7^ 0,/i8 7^ 0). The magnetization of the superconducting CFL phase is studied as a 



function of eB / ji^ for ji 
are observed. 



500 MeV in 



23], where the same vAdH oscillations as appears in Fig. 4 
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FIG. 4: The dependence of the product M • B on the magnetic field of fixed chemical potential a) /i = 250 
MeV in the xSB phase and h) (x — 460 MeV in the CSC phase. The linear regime in both phases starts at 
SB -0.45- 0.50 GeV^. 



In what follows, we will first study the ^-dependence of as and A^. We then present the phase 
diagram /ic — eB of the 2SC quark matter at zero temperature. Let us start with the case of zero 
magnetic field. In Fig. 5, the /^-dependence of (Tq in the xSB phase, as well as Aq and //g in the CSC 
phase are plotted for zero magnetic field. For our specific choice of free parameters k,Gs and Gd, 
(To = 323.8 MeV for ^ < ^c- Here, the critical chemical potential ^c = 325 MeV and the value of Aq 
for u ~ ^c is An = 78.0 MeV. Our results coincides qualitatively with the numerical results presented 
in [a] (see also 39| for a recent investigation of Cooper-pairing in NJL-type models). ^^ 
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In [g|, the quark mass mo ^ and therefore a mixed phase appears for /i > 340 MeV. 
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Analytical parameters 


Numerical fit parameters 


1] in % 


iB (GcV2) 


a (MeV2) 


h 


c (MeV^) 


a (MeV2) 


h 


c (MeV2) 


ria 


% 


??c 





8.49 X lO'^ 


0.19 2.29 X 10^ 


7.84 X 10^ 


0.13 


2.52 X 10^ 


8 


38 


10 



TABLE II: Numerical fit data for AQ(/i) from (|V.4p . The numerical values of the parameters arising from our fit 
are in good agreement with the expected analytical values arising from (IIV.39P [see rja, r/i, and rjc with ry defined 
in (rv3]) ]. 

FIG. 5: The /x-dependence of ctq in the xSB phase, and Aq in the CSC phase for eB = (left panel). The 
/i-dependence of /ig for eB = (right panel) . 



We can compare the /i-dependence of Aq arising from our numerical calculation with the relation 
PV.39P arising from our analytical results for vanishing magnetic field. To do this we have fitted our 
numerical data with a function 



^o(^) = {a-bfi )exp( ^ 



(V.4) 



similar to ()IV.39p . Here, a,b and c are free parameters. The numerical values of these parameters 
arising from our fit are in good agreement with the expected analytical values arising from pV.39p 
(see Fig. 6 and Table II). This can be quantified by defining 

Analytical value — Numerical value 



Yj = - (V.5) 

(Analytical value + Numerical value) /2 

as a measure for the variation of the numerical value with respect to the average of analytical and 
numerical values. In Table II, T/a)% and r]c are less than 50%. Note that the difference between the 
analytical and numerical values of a, b and c lies on the approximations that are made to determine 
analytically Ag in (jIV.SOp . 
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FIG. 6: The dots are the numerical values of A§. The solid line is the corresponding fit of Aq(/x) from (jV.4|) . 
The fit parameters a, b and c are listed in Table II. The regression parameter R^, as a measure of reliability of 
the numerical fit, is in this case R^ — 0.999852. 



Let us now concentrate on the case of non-vanishing magnetic field. In Table III, we have summa- 
rized our numerical results for critical chemical potential fic, the mass gap ctb for ^ < ^c and the 2SC 
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gap /S.B at /i ~ ^c- The critical chemical potential /Xc and the ^SB mass gaps (TB(/i < ^c) increase by 
increasing the external magnetic field. In the vicinity of the phase transition from xSB to CSC phase, 
the CSC mass gap AsifJ- — ^c) also increases by increasing the magnetic field. 



iB in GeV2 


/ic in MeV 


o'sip < ^J'c) in MeV 


AsifJ. ^ Mc) in MeV 


0.002 


316.0 


312.0 


70.0 


0.005 


287.0 


311.0 


61.0 


0.01 


324.2 


322.4 


78.0 


0.04 


322.3 


326.4 


77.6 


0.05 


321.9 


327.8 


77.9 


0.10 


311.3 


338.4 


76.9 


0.20 


296.0 


366.0 


50.8 


0.30 


315.0 


420.7 


73.2 


0.40 


329.0 


464.8 


80.4 


0.46 


329.6 


477.0 


82.6 


0.48 


334.0 


487.0 


85.7 


0.50 


340.8 


499.0 


88.8 


0.60 


392.4 


593.0 


116.9 


0.70 


450.7 


700.0 


147.7 



TABLE III: Numerical results for critical chemical potential ^c, the mass gap a-g for /i < /if, and the 2SC gap 
As at /Lt ~ lie- 

The ;U-dependence of as and A^ are presented also in Fig. 7. There is a first order phase transition 
from the xSB to the CSC phase [see also Fig. 9 for more detail on the phase structure in Hc — eB 
)lane]. Because of our specific choice mo = and Gd < Gs, no mixed broken phase appears at fj, > fic 
9|, and the xSB mass gap (Tb(//) is constant for /i < jic- For small value of eB, the CSC mass gap 
AsifJ-jcB), is increasing with /i. The magnetic field enhances the chiral symmetry breaking. This 
is known as the phenomenon of magnetic catalysis 31( , which is also observed in [35(] . In the linear 
regime, i.e. for eB > 0.45 GeV^, A^ is decreasing with /x. To study the linear regime in detail, we 
have fitted our numerical data for A^ as a function of /x and fixed eB, with a function similar to 



A 



{a-bf/), 



(V.6) 



where a and b are free parameters, that depend on eB. In Table IV, we have compared the expected 
analytical results for the parameters a and b, with the corresponding results from fitting our numerical 
data with ()V.6P for different eB. For eB > 0.44 GeV^, rja and rjh are less than 50%. 



In Fig. 8, the /x-dependence of the color chemical potential /xg for different eB is demonstrated. As 
in the previous cases, we expect that, in the linear regime eB > 0.45 GeV^, the fi dependence of /xs is 
given by a function similar to pV.30p . that arises analytically in the LLL approximation. We define 
therefore a function 



f^8 



a + biJ? 



(V.7) 
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FIG. 7: The /.t-dependence of as in the xSB phase (red hnes), and As in the CSC phase (blue hnes) for different 
values of eB. The xSB mass gap (T_b(/^, eB) is constant in /i < /ic and increases for increasing eB. The critical 
chemical potential fj,c increases for increasing eB. For our specific choice of parameters (ttiq = 0, Gd < Gs) 
no mixed phase appears. The CSC mass gap Ab exists therefore only at fj, > fic- The slopes of the curves 
appearing at /i > /ic are decreasing for increasing eB. The first order nature of the phase transition between 
xSB and CSC phases is visible. The dependence of the xSB and CSC gaps for vanishing magnetic field is also 
considered here to have a comparison with the /^-dependence of the gaps for non- vanishing eB. 

with arbitrary, ei?-dependent parameters a and b. In Table V, we have compared the data that arise 
numerically by fitting the numerical values of fig with (IV. 7p for different eB. As in the previous case, 
the difference between the numerical fit data and the expected analytical values of a and b arising 
from pV.30P minimizes in the linear regime for eB > 0.40 GeV^ {rja and 7]^ in Table V are less than 



Finally, we will present the phase structure of the model in a /ic — eB plane in Fig. 9. In particular, 
we are interested on the effect of the color chemical potential /tg on the phase structure of the model. 
In Fig. 9a (Fig. 9b) the phase structure for /ig = (/xg 7^ 0) is plotted. Because of our specific 
choice of parameters, we expect xSB and CSC phase without mixing. A normal phase can also exist. 
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Analytical parameters 


Numerical fit parameters 


77 in % 


eB (GeV2) 


a (MeV2) 


b 


a (MeV2) 


h 


Va 


rib 


0.04 


5.27 X 10-24 


1.32 X 10-28 


-1.06 X 10^ 


-0.077 


200 


200 


0.10 


1.62 X 10-6 


1.62 X 10-" 


-8.37 X 102 


-0.077 


200 


200 


0.30 


1.91 X 102 


6.37 X lO-'' 


-4.53 X lO-"^ 


-0.094 


218 


201 


0.40 


2.27 X 10^ 


0.006 


+4.17 X lO-"^ 


-0.018 


59 


400 


0.44 


4.53 X 10^ 


0.010 


+7.38 X lO-"^ 


+0.008 


48 


22 


0.46 


6.14 X 10^ 


0.013 


+8.47 X lO-"^ 


+0.015 


32 


14 


0.50 


1.05 X 10'' 


0.021 


+ 1.09 X 10* 


+0.025 


4 


17 


0.60 


3.03 X 10* 


0.050 


+2.26 X 10'' 


+0.058 


29 


15 



TABLE IV: Numerical fit data for A^ as a function of /i from (jV.6p . In the linear regime, i.e. for eB > 0.45 
GeV2, the numerical values of the parameters arising from our fit are in good agreement with the expected 
analytical values of the parameters from (|IV.33p [see rja and rjb with 77 defined in (|V.5[) ]. 
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FIG. 8: The /i-dependence of the color chemical potential /zg as a function. The numerical data for /^g are 
fitted in Figs 8d — 8/ by (IV. 71) . In the linear regime, the fitted curves [solid (blue) lines in 8d — 8/] are in good 
agreement with our numerical data. For eB — 0.44, 0.5 GeV2 and 0.6 GeV2, the regression parameter E? ^ as a 
measure of rehability of numerical fits are R^ = 0.999991, 0.999924, 0.999952, respectively 

where the mass gaps ub and A^ corresponding to xSB and CSC phases vanish identically.^^ To check 
this, we consider the gap equations and the color neutrality condition (|IV.ip . We have looked for the 



^^ As it is known from [5], in the regime of large chemical potential, jj, > 500 MeV, the 2SC phase goes over into the 
three-flavor CFL phase. In the present two-flavor model, we only assume that a normal phase may exist, and, if so a 
phase transition will occur from the color superconducting 2SC phase into this normal phase (see Fig. 9). Hence, the 
present results concerning the transition from CSC to the normal phase is only of theoretical nature. To include the 
CFL phase, we have to extend the model to three-flavor superconductivity including up, down and strange quarks. 
This is indeed beyond the scope of the present paper and is planned for future publications. 
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Analytical parameters 


Numerical fit parameters 


ri in % 


eB (GeV^) 


a (MeV2) 


h 


a (MeV2) 


h 


Va 


'/fe 


0.04 


1.80 X 10^ 


6 


-1.01 X 10^ 


+7.53 


207 


23 


0.10 


4.50 X 10^ 


6 


-4.93 X 10^ 


+ 161.14 


204 


186 


0.30 


1.35 X 10^ 


6 


+3.04 X 10^ 


+13.04 


77 


74 


0.40 


1.80 X 10^ 


6 


+ 1.93 X 10^ 


+6.27 


7 


4 


0.44 


1.98 X 10^ 


6 


+2.34 X 10^ 


+4.11 


17 


37 


0.46 


2.07 X 10^ 


6 


+2.43 X 10^ 


+3.89 


16 


43 


0.50 


2.25 X 10^ 


6 


+2.46 X 10^ 


+4.49 


9 


29 


0.60 


2.70 X 10^ 


6 


+2.45 X 10^ 


+5.35 


10 


11 



TABLE V: Numerical fit data for /ig as a function of ji from (IV. 71) . In the linear regime, i.e. for eB > 0.40 
GeV^, the numerical values of the parameters arising from our fit are in good agreement with the expected 
analytical values of the parameters from (IIV.30I) [see 'qa and rjb with r] defined in (jV.5[) ]. 



global minima of the system in two different regimes: ^c — 350 — 450 MeV and /ic — 750 — 800 MeV. 
As it turns out, in the first regime corresponding to ^c — 350 — 450 MeV, the minima of il^ft from 
pil.24p are given by {as / 0, A^ = 0, /^s = 0) for ^ < /^c as well as {as = 0, A^ / 0, /^s / 0) for 
IJ, > He- In the second regime corresponding to ^c — 750 — 800 MeV, however, the global minima are 
{ub = 0,Ab 7^ 0,/.i8 7^ 0) for /i < /Uc, as well as {as = 0, A^ = 0, /is = 0) for ^ > ^c- We conclude 
therefore that a phase transition from xSB to CSC phase occurs in the first regime at ^c — 350 — 450 
MeV, and a phase transition from CSC to the normal phase occurs in the second regime ^c — 750 — 800 
MeV. In the following, we denote the value of O^g at the global minima by il^sB = ^off('7B, 0, 0; ^u, eB), 
r^csc = f^cft(0, Ab, fJ.8', fJ-, eB), and ilNo^nj^i = ilcff (0, 0, 0; fi, eB) corresponding to the xSB, CSC and the 
normal phase, respectively. For different values of {fj,,eB), the xSB phase is defined by fi^sB < ^csc 
and the CSC phase by Jlcsc ^ ^xsb- Moreover, the exact value of //c for the first order phase transition 
from xSB to the CSC phase [the lower (red) solid line in Fig. 9a and 9b] and from the CSC phase 
to the xSB phase [(green) solid line in Fig. 9b] are then defined by 0,^s^ = r^csc and Clcsc = ^Normal, 
respectively 41[. As for the second order phase transition between the CSC and the normal phase, an 
analysis similar to [TJ] is performed. In Fig. 9a, the phase diagram of the 2SC model in a Hc — eB plane 
is plotted for vanishing /ig. A first order phase transition occurs between the xSB and the CSC phase 
in the regime fic — 350 — 450 MeV (solid red line). This confirms the results by 201, where a first order 
phase transition is observed for fixed value of eB = 0.05 GeV^, and various Gd/Gs- The transition 
from the CSC into the normal phase is of second order and occurs in the regime fj,c — 750 — 800 MeV 
(dashed black line). According to the phase diagram for nonvanishing /xg in Fig. 9b, however, whereas 
the transition from the xSB to CSC phase is still of first order (solid red line), the second order phase 
transition for small values of eB (dashed black line) goes over into a first order phase transition at 
fic ~ 755 MeV and eB ~ 0.13 GeV^ (solid green Hne). Moreover, at eB ~ 0.4 GeV^, ^c suddenly 
decreases and increases once again by increasing the external magnetic fields in the CSC regime. The 
CSC regime is nevertheless suppressed in the linear regime eB > 0.45 GeV^ by the external magnetic 
field (see Fig. 9b). 
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FIG. 9: The phase diagram of 2SC model is plotted in a /Xc — eB plane for a) fig = and b) fig ^ 0. The solid 
(red) lines in a) and b) indicate first order phase transitions between the xSB and the CSC phase. The dashed 
(black) lines in a) and b) are second order critical lines between the CSC and the normal phase. As it is shown 
in b), for /is 7^ and at ^c — 755 MeV and eB ~ 0.13 GeV^, the second order phase transition goes over into a 
first order phase transition between the CSC and the normal phase [solid (green) line]. At eB ~ 0.4 GeV^, Hc 
suddenly decreases and increases once again by increasing the external magnetic fields in the CSC regime. 

VI. CONCLUDING REMARKS 

In this paper, we have studied the effect of constant magnetic fields on the formation of bound states in 
the chiral as well as the color symmetry breaking phase. In the first part of the paper, after introducing 
a two-flavor NJL type model including meson and diquark condensates, we have computed the one- 
loop effective action and the thermodynamic potential of the theory at zero temperature and finite 
density. Neglecting the quark mass mo and choosing the diquark-to-chiral coupling ratio Gd/Gs < 1 
[9|, we can consider the xSB and CSC phases separately. The xSB and CSC mass gaps aB and A^ as 
well as the color chemical potential /xg are determined analytically in the limit of very strong magnetic 
fields. In this limit, the dynamics of the system is dominated by the lowest Landau level and therefore 
the effect of all higher Landau levels are negligible. According to [3ll], in this limit, as a result of 
dimensional reduction from D to D — 2 dimensions, the formation of bound states and consequently 
a dynamical symmetry breaking will be possible even for weak interactions between two fermions. 
This is the phenomenon of magnetic catalysis, discussed widely in the literature in the past few years 
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34( 1 . Denoting the dimensionless coupling constants in the xSB and CSC phases by gs ~ G^A^ 
and gd ~ G^A'^, we have determined the mass gaps for different regimes of gg and gd- Here, A is 
certain momentum cutoff. In [g], the xSB and CSC mass gaps of a similar 2SC model was determined 
numerically for vanishing magnetic field. Introducing a large momentum cutoff A and performing 
appropriate approximations, we have determined analytically the mass gaps ctq and Aq as well as /xg 
corresponding to xSB and CSC phases for zero magnetic field too. 

In the second part of the paper, a detailed numerical analysis is performed to explore the effect of 
any arbitrary magnetic field on the mass gaps a and A and the color chemical potential /ig. First, the 
dependence of a and A on various eB G [0, 1] GeV^ is plotted for fixed fi = 250 MeV and fi = 460 
MeV in the xSB and CSC phases, respectively. For small values of eB, we observe small van Alfven-de 
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Haas oscillations, that appear, according to [221 . |23[ also in the CFL phase for /i = 500 MeV. Same 



oscillations appears also in the magnetization M for the same fixed value of chemical potentials. At 
eB ~ 0.4 — 0.5 GeV^, the oscillations end up in a linear regime, where we believe that the dynamics 
of the system is described exclusively by the LLL. This can be checked by comparing qualitatively 
the numerical dependence of as and A^ for eB > 0.45 GeV^ and fixed /j,. The //-dependence of a 
and A are then considered for various eB. Our numerical results for vanishing magnetic field coincide 
with the numerical results presented in [6J]. The ^-dependence of as and As as well as ns are then 
considered for various finite eB. The numerical results in the linear regime, i.e. for eB > 0.45 GeV^ 
are comparable with our before mentioned analytical results in the limit of large eB. This is shown 
using appropriate numerical fits. The phase structure of the xSB and CSC phases in a, fic — eB plane 
is also presented. We are in particular interested on the effect of the color chemical potential //g on 
the phase diagram of the model. For both ^g = as well as /ig ^ 0, the transition from the xSB phase 
into the CSC phase is of first order, and occurs in the regime Hc ~ 350 — 450 MeV and for eB G [0, 0.7] 



GeV^ (see Fig. 9a). This confirms the result in 26|], where a first order phase transition is observed 
between the xSB and CSC phase for fixed eB = 0.05 GeV^ and various Gd/Gs ratios. Assuming 
that the CSC phase goes over into a normal phase at ^ > 500 MeV,^^ it turns out that whereas for 
//g = 0, this transition is of second order, for nonvanishing /xg, a second order phase transition occurs 
first for small eB. It goes then over into a first order phase transition at /ic — 755 MeV and eB ~ 0.13 
GeV^. At eB ~ 0.4 GeV^, fic suddenly decreases and increases once again by increasing the external 
magnetic fields. The CSC phase is nevertheless suppressed in the linear regime eB > 0.45 GeV^ by 
the external magnetic field (see Fig. 9b). 

At the end, let us just emphasize that the study of color superconductivity in the presence of 
constant magnetic fields has not only astrophysical consequences in forming the structure of compact 



star cores, it may be also relevant for future heavy ion collision experiments. Recently in 17|, the 
accessibility of the color superconducting 2SC phase in the heavy ion collisions is investigated. It is 
stated that for high enough collision energies the 2SC may be accessible in future collision experiments. 
On the other hand, there are various evidences of the creation of very strong magnetic fields in non- 
central heavy ion collisions [15|, ll6|. It would be interesting to perform similar analysis as in 17] 
considering the presence of constant magnetic fields. To do this, the effect of finite temperature on 
the phase diagram of the 2SC superconducting phase in the presence of constant magnetic fields has 



also to be considered. This will be reported in future publications [111 ]. 
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